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FOREWORD 

This  report  is  submitted  as  part  of  Contract  AF 
33  (038)  21406,  whicn  was  initiated  by  the  Office  of  Air 
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School  of  Aeronautical  Engineering  at  Cornell  University 
in  March  1951.  The  report  is  one  of  a  series  to  be  published 
as  the  result  of  the  work  carried  out  under  this  contract. 

The  present  investigation  was  carried  out  partly 
under  the  support  of  the  contract  mentioned  above  and 
partly  as  the  author’s  thesis  investigation  for  the  degree 
Master  of  Aeronautical  Engineering.  The  subject  of  the 
investigation  was  suggested  to  the  author  by  Professor  T 

J.  M.  Wild,  now  Acting  Chief  Engineer  of  the  ARO  Corporation. 

The  author  wishes  to  acknowledge  his  indebtedness  to  Mr. 

Wild  and  to  the  Faculty  of  the  Graduate  School  of  Aeronautical 
Engineering  of  Cornell  University  and  in  particular  to 
Professor  W.  R.  Sears  for  his  suggestion  of  this  study 
and  his  invaluable  advice  and  encouragement  throughout 
the  study  and  to  Professor  N.  Rott  for  his  most  helpful 
suggestions  and  criticism  while  preparing  this  report. 

The  work  on  this  project  was  conducted  under 
Research  and  Development  Order  No.  465-5-6,  Practical  Problems 
in  Aerodynamics.  Mr.  L.  S.  Wasserman  of  the  Plight  Research 
Laboratory,  Wright  Air  Development  Center,  was  the  project  engineer 
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ABSTRACT 


The  flow  characteristics  behind  weak  stationary 
shock  waves  reflected  from  various  edges  lying  in  the  main¬ 
stream  direction  are  determined  according  to  linearized 
theory.  Five  different  types  of  edges  axe  considered, 
made  up  of  various  combinations  of  solid  and  free  plane 
surfaces.  An  assumption  regarding  the  singularity  of  one 
of  the  perturbation  velocity  is  required  in  order  to  render 
the  solutions  unique. 

By  superposition  of  such  basic  "edge"  solutions, 
the  flow  behind  a  snock  wave  reflected  from  a  wall  with  a 
slot  (i.e.,  a  strip  of  free-surface  between  panels  of  solid 

m  * 

wall)  and  from  multiply- slotted  walls  are  obtained.  These 
*  solutions  apply  only  to  regions  upstream  of  multiple  inter¬ 

actions  of  the  slot  edges;  however,  these  regions  include 
the  most  interesting  regions  of  flow  in  the  case  of  re¬ 
flection  from  slotted  wind-tunnel  wall,  for  example. 

The  relation  of  the  single- slot  problem  to  the 

problem  of  a  narrow  rectangular  supersonic  wing  is  discussed. 

PUBLICATION  REVIEW 

The  publication  of  this  report  does  not  constitute 
approval  by  the  Air  Force  of  the  findings  or  the  conclusions 
,  contained  therein.  It  is  published  only  for  the  exchange  and 

stimulation  of  ideas. 
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INTRODUCTION 


In  recent  supersonic  wind  tunnel  research  work,  there 
occurs  sometimes  the  necessity  of  a  slot,  or  multiple 
slots,  along  the  test  section,  parallel  to  the  stream. 

The  slot  is  usually  in  the  form  of  a  strip  of  free  surface, 
say,  stationary  air,  between  neighbouring  solid  walls. 

Also  there  may  be  mutually  perpendicular  edges  in  various 
combinations  of  solid  wall  and  free  surface.  When  a  two- 
dimensional  shock  wave  from  a  wind-tunnel  model  is  re¬ 
flected  from  the  slot  or  the  edge,  the  flow  character¬ 
istics  after  reflection  becomes  three-dimensional.  The 
determination  of  the  flow  field  is  the  subject  of  the 
present  study. 

In  Part  I  of  this  thesis,  basic  ideas  and  equations 
are  discussed  and  derived.  In  Part  II,  flow  character¬ 
istics  after  a  weak  shock  hitting  an  edge  are  determined. 
(Five  different  kinds  of  edge  problems  are  considered). 

Then  the  problems  of  a  shock  reflected  from  a  single  slot 
in  infinite  solid  wall  and  from  a  multiply-slotted  solid 
wall  are  considered  in  Part  III  and  IV  respectively. 

In  Part  V,  correspondence  of  boundary  conditions  for 
velocity  component  w  of  the  narrow  rectangular  supersonic 
wing  and  u  of  the  narrow  slot  problem  are  established. 

It  is  shown  then  that  their  difference  in  singularities 
makes  the  solutions  of  the  two  problems  not  identical. 
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PART  I 


BASIC  CONSIDERATIONS 


Let  U  be  the  undisturbed  supersonic  velocity  along 
x  axis.  When  the  stream  passes  over  some  boundaries  or 
obstacles,  there  is  velocity  disturbance  which  may  be 
denoted  by  its  components  u,  v,  w  along  x,  y,  z  axes 
respectively.  Assume  the  flow  to  be  irrotational  and 
steady  and  denote  <J>  as  the  perturbation  velocity  potential 
such  that 


U-2t_ 

■9X 


v  _ 


w 


IE 


These  perturbation  velocity  components  as  well  as  pertur¬ 
bation  density,  pressure  are  supposed  so  small,  compared 
with  the  undisturbed  values,  that  their  square  terms  or 
cross  product  can  be  neglected.  Using  Euler's  equations 
of  motion,  continuity  equation,  together  with  the  equa¬ 
tion  of  state  of  gas,  assumed  perfect,  one  can  easily 
derive  the  well-known  Prandtl-Glauert  equation 


■****  -^31  -  na»  "*  °  (le 

Where  xx.  denotes  q>  ,  u,  v,  w  or  perturbation  pressure, 
and  m^  =  -1,  M,*  being  the  free  stream  Mach  number. 

Before  going  into  the  edge  problems,  it  will  be 
advisable  to  present  first  the  two  dimensional  result  of 
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a  shock  reflected  from  a  solid  surface  or  from  a  free 
surface. 

Assume  a  weak  compression  shock,  such  as  that 
results  from  a  supersonic  flow  about  a  thin  airfoil 
at  a  small  angle  of  attack  c<,  impinges  on  a  solid  wall, 
infinite  in  extent0  A  compression  shock  will  be  re¬ 
flected  from  the  wall.  Consistent  with  the  assumption 
of  small  perturbation,  the  angle  of  reflection  will  be 
same  as  angle  of  incidence  which  should  be  taken  to  be 

equal  to  Mach  angle  sin”^  —i  ol*  tan"^  — i  in  the  present 

M «»  m 

approximation,.  Values  of  perturbed  velocity  components 
are  summarized  in  the  following  figure: 


If  a  weak  compression  shock  impinges  on  a  free 
surface,  an  expansion  wave  will  result  after  reflection. 
Perturbed  velocity  components  has  the  following  values: 
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Again,  within  the  present  approximation,  both  angle 
of  incidence  and  angle  of  reflection  will  be  equal  to 
Mach  Angle.  Moreover,  although  the  free  surface  will 
deflect  downward  through  an  angle  of  after  reflection, 
we  will  still  consider  it  undeflected  in  the  edge  problems 
discussed. 

Five  different  kinds  of  edges  are  going  to  be  con¬ 
sidered  in  Part  II: 


(A) 


- { -  77777777777 

free  surface  solid  wall 


(1.4) 


(B) 


(C) 


(D) 


(E) 


'97777777777777777 77 
v  *  solid  wall 

^solid  wall 
4^ 


1 777777777777777777 
|  *  solid  wall 

I 

I  free  surface 
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^  ^free  surface 

/ 

✓ 
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I 

!  free 
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£*free  surface 
surface 


(1.5) 


(1.6) 


(1.7) 


(1.8) 
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When  a  shock  hits  an  edge,  the  flow  problem  around 
the  edge  becomes  a  three  dimensional  one.  However,  due 
to  the  principle  of  forbidden  signals  of  supersonic  aero¬ 
dynamics,  (Reference  1),  the  edge  effect  is  restricted 
within  the  Mach  cone  which  has  an  apex  at  the  intersection 
point  of  the  shock  and  the  edge.  Outside  Mach  cone,  the 
flow  characteristics  are  still  two  dimensional  as  mentioned 
before.  Furthermore,  the  present  problems  contain  no 
characteristic  length.  By  dimensional  theory  argument 
(Reference  2),  one  concludes  that  the  edge  effect  results 
in  a  conical  problem,  i0e.  all  flow  characteristics  are 
constant  along  the  ray  radiating  from  the  apex  of  Mach 
cone.  The  powerful  technique  of  conical  flow  may  then 
be  employed  here. 

Change  cartesian  coordinates  x,  y,  z  to  cylindrical 
one  x ,  vs'  ,  ix>  where 


Or  =  J  ^  -+  ? 

to  *=  ton-*  -fr 


(1.9) 


Equation  (1,1)  reduces  to 

frl  ~  tx/zcy  —  —  .O  .>*.  —  JL 

WWT  ■x*-YTy  m-v 


t Tr  -"■"'■ay  ro-v  =»  o 


JX. 


(1.10) 


If  denotes  the  cartesian  velocity  components 
u,  v,  w  or  perturbation  pressure,  but  -not  velocity  potential .  <t> 


„G 

■<1: 
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conical  flow  properties  ensure  that  XI  will  be 
of  ^  ,  M  only.  By  introducing 

tj  _ 

1  x 

equ0  (1.10)  becomes 

(mY-  0^09  +  (2m'y-±.)&?  -±awc0  m  o  (lel2> 

Use  Tschaplygin  transformation  (Reference  2) 

<r  n  1  3 

^jTT^r  U.i3) 

(Note  that  on  Mach  cone:  fj  =  1/m,  s  =  1) 

(1.12)  reduces  to 

•a35+^i+  ^0^=0  (1.14) 

which  is  the  Laplace  equation  of  SX  in  s,<o 

To  find  n.  -  u,  v,  w  within  the  Mach  cone,  one  has 
to  solve  equation  (1.14)  with  the  appropriate  boundary 
values  of  n  or  (the  normal  derivative)  described 
completely  along  the  boundary  (some  of  these  boundary 
values  are  immediately  known,  some  would  be  determined 
with  aids  of  irrotationality) .  Thus  we  usually  have  the 
so  called  "mixed  boundary  value  problems."  Technique 
to  find  the  solution  and  uniqueness  of  solution  and 
singularity  behaviors  will  be  considered  in  Part  II. 

Note  that  if  one  of  the  velocity  components  is  known, 


a  function 

(1.11) 
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the  other  two  may  be  found  by  irrotationality  condition, 
ecg.  if  u  is  known,  v,w  may  be  obtained  from 


=jt- 


co  cm. <o  *34 


)x  ”50) 


7  *7 


J  <^7  +  F(co)  (1.15) 


=  -J[ 


^  (1.16) 


or  in  (s,  id) 


V  o  in 
x 


jr^ 


3U  _  ivs 
•^co 


co  f(u>)  (1o17) 

J 


C-*S  +  (la8) 


Usually  one  is  interested  in  the  pressure  distribu- 
tionc  As  is  well-known  in  the  perturbation  theory,  the 
pressure  coefficient  Cp  is  given  by 


S  = 


t  -  K  _2  JL 

±UK  *  * 


(1.19) 


WADC  TR  52-161*. 


6 


PART  II 


SHOCK  REFLECTED  FROM  AN  EDGE 

We  shall  consider  the  flow  characteristics  of  the 
edge  problems  for  the  five  cases  mentioned  on  page 
one  by  one„  In  case  At  we  shall  develope  the  method 
solving  the  problem  in  detail,  and  find  the  three  velocity 
components  u,  v,  W  completely,.  In  other  cases,  only  u 
(i„e0  the  pressure)  is  solved;  v  and  W  can  always  be 
found  by  means  of  (1.15)  —  (1 0 18 )  if  required. 

The  notation 

K  5  -  W 

m 

or  mK  =  -  Uot 
is  used  in  this  part. 

Case  A: 


Ac  is  2K  and  along  BC  and  Bo,  zero.  Along  oA  ,  we  only 
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have  W=oc  However,  as  there  is  no  change  of  W  along  x 
direction,  we  have  =  o  and  hence  by  irrotationality 
*  o  along  OA0  Thus  the  problem  of  finding 
u  due  to  edge  effect  is  to  solve 

rr*  ,,  _  i  i 

V_  U  s  . .  ■+  —  —  ■  + - —  ■*«  o 

■  s.*o  /&3*  5  "as  '%*>*■  (A02) 

in  the  semi-circular  region  subjected  to  boundary  condi¬ 
tions  as  shown  in  .Figc  ( A0 3 ) : 


To  solve  (Ac3),  use  conformal  transformation 

i. 

5  =  «k*- 

where  _  _  ie 


s  -  r 


r  =  s 


The  problem  is  transformed  to 


( A  o  4 ) 

(A05) 

(A.6) 

(A. 7 ) 
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As  the  sides  with  u  =  o  and  —  =  o  requires  antisym¬ 
metric  and  symmetric  continuation  respectively,  Fig„  ( A0 7 ) 
may  be  completed  to  a  full  circle: 


(A08) 


One  of  the  fundamental  solution  of  Equation  (a.2) 
is  u=£.,p  ,  where  is  the  angle 
shown  in  the  Figure.  This  may 
be  verified  by  direct  substitu¬ 
tion  and  may  be  imagined  as  the 
potential  of  a  vortex  of  strength 
-k  situated  at  ,  Q„  „  To  solve 


may  place  vortices  of  strength  + 


x 


xtc 

X 


and 


at  A.,  a,  c,  D  respectively.  Then  at  any  point 


r  ©  inside  the  circle, 


7T  a. 


TT 


+*  & 


IF  ®4- 


constant 


-  f  (•»-©»)  - +  constant 
— .  *^-(  Ai  ~  A*)  +  constant 

To  determine  the  arbitrary  constant,  take  point  pc  Here 
and  u=o  as  givenc  Hence  the  constant  has 

to  vanish. 
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It  can  be  easily  seen  that  the  boundary  conditions 
along  AB,  BC,  CD  are  all  satisfied.  For  example,  along 
AB,  A9  arid  A  2  always  take  constant  values  ^  and  -^Ir¬ 
respectively.  Then  u  =  2K  as  required. 

To  express  A^,  A  2  ^eITns  of  r,  ©  ,  remember 
that  the  circle  is  an  unit  circlec  By  some  algebric 
manipulation,  we  have 


fan  A,  - 

2.  sln<f  (Cos  —  v-  c<y*  ©) 

Cos  a»j>  — i  rcos  cos©  +  ir1 

( A.10) 

-ten  Ax  - 

Z  aln  ^  (  Cos  -+-  y  Cos  &  ) 

(A. 11) 

Cosaf  4-  arcos  ©  cos^  4-  r3- 

A  s  A,  —  Ax 

_  1  -I  C  4  Sih  t  —  ra)  Cos  8 

(  I r4  +  (  Cos  -2  cos*©)  -v  | 


f 


(A012) 


Here  ^  =  ~ 


A  -  tsrT1  C  l"tl~  1 

<  l"*  —  4-  t'1Cos1*  I  J  ( A,  13 ) 

By  means  of  (A, 6) 

a  4.  f  x  £3  <■  t-  s)  7 

A  “  ]  S*  -45  Co^-^-  +  |  j  ( A,  14 ) 

Finally,  by  equ.  (1.13)  the  solution  may  be  expressed 
in  physical  coordinates  b  ,  «*: 
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where 


( Ac  15 ) 


a- got,  A_ 
«n  '  7^ 


A.  =  t"on 


^  f  1-^9)  Co: 


co 


1  —  xmr) 


Co* 


and 


A  < 

7T 

wKen 

a  i  to  *  7t  , 

A  = 

X 

when 

0  &  to  &  ^  t 

tt0- l 

A.  n 

0 

when 

2 

m  *)*  I 

(A, 16 ) 


(A»17 ) 


The  pressure  coefficient  is,  due  to  (1.19) 

C„  _  _  _A_ 

”  **  *  TV  ( A « 18 ) 

where  A  takes  same  value  as  in  (A. 16),  (A„17)„ 

On  free  surface:  to  =  -7c  3  a  =  o 

o 

on  solid  surface:  co  »  o  ,  K  =»  Co»“l  (  1-  am  ^  J 

CP  “  ^  cos-'O  _  iWv,  )  (a, 19) 

The  distribution  of  cp  along  these  surface  and  a  set  of 
constant  pressure  lines  (isc-bsr)  has  been  plotted  in 
Figure  1  and  2  (P0  47*  46)  respectively,, 

The  point  c  in  Figure  2,  p0  4g  is  quite  interesting. 
At  this  point  the  boundary  values  given  by  the  two- 
dimensional  results  are  discontinuous,  and  the  expression 
of  (Ac.16)  takes  the  indeterminate  form  tan”^  (~).  It 
follows  that  every  iso-bar  of  different  values  of  cp 
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converges  to  the  point.  Moreover,  the  tangents  of  all 
iso-bars  are  horizontal  (parallel  to  boundary  surface) 
at  c  except  the  one  of  Cp  =  2 -7^  ,  which  makes  135°  with 
the  horizontal  (see  Appendix  I), 

The  position  variable  used  in  Figure  2  is  7 „ 

However,  if  we  assume  m  =  1  (i„e0  M&>  =  J~2 )  and  give  x  a 

definite  value,  say  1,  represents  the  actual  physical 
radiel  distance  It  is  interesting  to  note  that:  when 

going  inward  radially  from  points  along  AC  the  pressure 
decreases  quite  rapidly,  and  from  points  along  BC  the 
pressure  increases  very  rapidly.  The  phenomena  becomes 
more  acute  in  the  neighborhood  of  Ca 

These  statements  may  be  clearly  seen  from  the  following 
expression  of 


M. 

'bW  7\ 


(  l  -  A,m>)  /W|) 

^ 0  (  i  -  0  )  (.  l  —  nT  rj1-  coj 


(A, 20) 


Here,  we  find 


and  when  >n»)  =  1 


if  l 


— ^  <  0  if 


>  0  i.e.  4(H^<7==r 

J  Jim  n 

;  i.e.  _1_ 


_ 


-for  o-s,  to 

2l 


2Sf_  „ 

■"6  to-' 


-for  -^  <  to  <  7T 
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From  these  considerations,  we  are  easily  led  to  assert 
that  AC  and  BC  are  expansion  wave  and  compression  wave 
respectively,  They  may  be  considered  as  the  continuatic 
of  the  originally  reflected 
waves,  as  shown  in  the  accom 
panied  figure. 

Now  let  us  consider  V  0  On  Figure  (A.l),  along 

AC  and  BC,  V  ~  o  „  Along  oA,  =  o  gives  -  —  ,  o  _ 

Alonj  BO,  —  o  .  Hence,  -^*  =  o  b_u  irro+a-Won&li'f 

i,e„  V  is  constant  along  x.  Since  v=  o  on  Mach  cone 
(bounded  by  2-dim,  flow),  this  constant  vanishes,  and 
we  have  v =  o  on  the  entire  free  surface.  The  problem 

of  solving  V  is  to  solve  =  o  in  the  following  region. 
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By  (A.15) ,  (A. 16) 


-au  ^  _  3.Tfd>  Oxi  -x  (  i  —  a-*1?  '^,3~ jgr) 

'  tt  J  n»*7  ( i-  (  t-  ^ 

til  ixcrt*  /al*  ~  (.  \+  z-'nn  )J  »>*'7(  i->»h 

ZU3  "'X  (  franco) 

Due  to  (1.15) 


V  =  ^.CoSJ5 


r  d(^o> 

J  (.»*»?)*  c i- 


n)i  (j- »>*■)) 


Hr  P(  CO  ) 


-  -  +  F(») 


When  **r)  =  1,  V  =  o  ,  Hence  F(W)  vanishes,. 


(Ac  22) 


(A. 23) 


_  AgsL  Co^  4-  -  22 


( A  o  24  ) 


Thus  V  is  not  identically  equal  to  zero,  but  has  a 
singularity  of  j=U  type.  This  solution  satisfies  all 
boundary  conditions  of  Fig.  (a021). 

For  W,  refer  to  Fig.  (A.l),  we  have  W=o  along  oA, 
AC  and  W  =  -2  trot  =  2  Km  along  CB.  Along  Bo,  since 
u  =  o  ,  v  =  o,  we  have  ^  =  o,  ||  =  o. 

By  means  of  continuity  equation  in  small  perturbation 


theory 


3*  1  'IVj 


7)W 

+  =o 


**  ^  ^  (A.  25) 

We  have  OX  =  =  o  along  Bo.  Thus  to  find  W  is  to  solve 

o  xr  'dh 
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(A. 26) 


The  problem  is  similar  to  Fig.  (A. 3),  to  which  it  can 
be  transformed.  Replacing  2K  in  (A. 9)  by  2Km  and  ui  in 
(A. 16)  by  (7r-ttj),  we  have  the  solution  for  problem 
Figc  (A.26): 

•  W  =  2 Km  -  2lXc*  —  (A.27) 

where  (  tlE.  (A. 28) 

(  -  4-S  -V  |  J 

=  w  (  ±IE5HZE5I ihM-l  (A. 29) 

(  ^  --sm1-—  £ 


and 


0  <  A  <  7V 

when 

0  Sr  U)  £  7T  , 

mq<.  \ 

e  « 

> 

1! 

when 

£  *  W  &  A  , 

I 

A  =  0 

when 

^  > 

mg  =  l 

There  may  be  also  some  singularities  which  satisfies 
Zelro boundary  conditions  (Note  that  Laplace  eqn0  is  a 
linear  differential  equation  and  superposition  of 
solutions  may  be  used).  Compute  W  by  equ.  (1.16). 
With  aids  of  (A. 22),  (A. 23)  we  have 


WADC  T ft  52-16U 


15 


w  =  -sin 


‘  +  ,(w) 
<  ^0  )*  t  <~  ">  9  J  *•  (  I-  *>y^«>\o) 

_  _  vlfg(  ^-l  f  •a-J  »)Q  (  I-  *t?t7~)  3^~  7 

*  |  i  -  ■s,h*-^  | 


(  A  .-27a 


-  —  J  -»?  +  *<“> 

-1 


Where  G(«>)  =  o  if  same  value  of  tan  (  )  QS  that 

in  (A. 30)  is  used.  (integration  detail  is  given  in 
Appendix  II) . 

Here  we  see  that  in  addition  to  the  solution  (A.27), 
we  have  a  solution  of  singularity  of  -=J=  type.  The  latter 
satisfies  the  following  conditions: 


(A. 31) 


By  direct  differentiation  of  v  and  W  (A.27aand  A. 24), 

it  can  be  easily  verified  that 

_  -aw  _ 

T>2:  'zy  ° 

is  also  satisfied. 


It  is  clear  now  that  to  solve  the  edge  problem  (and 


in  fact  many  other  problems  in  conical  flow)  one  should 
be  rather  careful  about  the  possible  existence  of  a 


singularity*  It  may  be  possible  to  set  up  the  boundary 
value  problems  separately  for  u,  v  and  w,  However,  the 
possible  existence  of  singularity  makes  the  solution  of 
Laplace  equation  not  unique  even  if  the  differential 
equation  and  the  boundary  conditions  are  properly  satis¬ 
fied.  To  make  solution  unique,  it  is  necessary  and  suf¬ 
ficient  to  make  physical  assumptions  on  singularities  for 
one  of  the  quantities  u,  v,  wc  In  our  case,  we  shall 
assume  there  is  no  singularity  in  u  (the  assumption  may 
be  justified  from  a  physical  point  of  view)  and  find  u 
uniquely  by  solving  the  mixed  boundary-value  problem, 
v  and  W  admit  singularities.  They  will  be  determined 
uniquely  by  means  of  irrotationality  conditions,  i„e., 
by  equation  (1.15),  (1.16),  if  u  is  known. 

It  is  also  interesting  to  notice  that  in  the  present 
case,  the  regular  part  of  W  is  a  mirror  image  of  u  (com¬ 
pare  A.iV^with  A. 15)  and  the  singular  part  is  a  mirror 
image  of  v  (compare  A.£.?awith  A. 24).  That  is, 


-  W (<),«0  -  m  u  (0,  -n-w)  +  vi>)  ,-n-cu) 


fi 


Case  B; 


fe+lec+ed  SkocK  wo vC. 


wahc  tr 


(B.l) 


17 


Owing  to  irrotationality  conditions  and  2-dimensional 


values,  the  problem  of  finding  u  becomes  that  of  solving 
the  following  boundary-value  problems: 


(B.2) 


Use  conformal  transformation, 


which  has  same  solution  as  that  of 
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To  solve  (B.5),  place  vortices.  K  and  -  K  at  A  and  B 

respectively.  For  any  point  (r  $  )  inside  the  circle, 

U  — ■  +•  constant 

as  JS +  constant 

A 


At  point  C#  A  *  9r*and  u  =  2K  as  given.  Hence 

* 

constant  should  be  equal  to  vYe  ^ave 

u - —  (  A  +  ^7r) 

To  express  A  in  r,  e^put  £  in  equ.  (A. 10) 


A  —  +an 


f'  ^  J3  (  >  -  *  F  Cos  »  ) 

)  ir1  —  ar  cose  -  | 


where 


?-  <  a  <  r  t 


when  o  £•  ftts 


C  o  <  A  <  7T 
Transform  back  to(s,  w,) 


i/vhe  n 


f  4  6  i 


(B.6) 


( Bo  7 ) 


(B.8) 


where 


JS(l 

-a5*c»S«.)  ] 

(B„9) 

Us* 

t  ,  f 

—  aS  CoS5«)  —  |  V 

?  < a  <  i* 

when  6  4  iO<  Y 

(B.10) 

e>  <  A  <  X 

when  2L  u)  ^  —  x 

-X  x 

The  pressure  coefficient 


(  A  +  -r  -re  ) 


(B.ll) 
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The  value  of  'ftp  along  the  solid  wall  has  been 
plotted  in  terms  of  physical  coordinates  ^  ,  w,  (See 
Fig.  3,  p.  49).  A  set  of  constant  pressure  lines  (iso¬ 
bar)  is  also  plotted  (See  Fig.  4.  p.  30). 

Here,  similar  to  what  we  have  in  case  A,  the  line  AC 
behaves  like  an  expansion  wave  and  BC,  like  a  compression 


Here  to  find  u,  one  has  to  solve 

U  =  ©  (C.2) 

with  boundary  conditions  given  in  Fig.  (C.3) 
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where  3  »re9,  r  =  se“  „  Problem  (C„3) 


which  has  same  solution  as  that  of 


Problem  (C„6)  may  be  considered  as  a 
of  the  following  two  problems:  (C,7)  and 


(C.3) 


(C.4) 

is  transformed 

(C.5) 


(C.6) 

superposition 
(C„8)  o 
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(C.8) 


These  two  problems  are  same  as  that  of  Fig,  (A. 8)  except 
herei}'-^  f-  for  (C.7)  and  (C„8)  respectively  and  boundary 
value  of  u  is  reduced  by  one-half.  Hence  by  (A„12), 
solution  of  (C.7)  will  be 


u  -  JL  W  (  - ? 

K  l  +  rxCi  -4  Cm*©)  +  |  5 


(C.9) 


where 


A 

<  IK 

when 

osotf 

r  < 

A 

— 

when 

£> 

r  « 

A 

—  0 

when 

cS 

r  * 

Solution  of  (C.8)  is 


x  r  cos© 
I  - 


) 


(ColO) 


(c.n) 


Where  o  s  ^  5  a.  ,  when  o  <  &  5 

So  back  to  problem  (C„3),  The  solution  in  (s,a>)  is 

y  _  XTo<  ~i^i 

^  -r r 


(C.12) 


where 


f  as‘(H  Cos  ^ 

V,  =*  +Qn  \  — x - 5“* - 1 - 

I  I  3*  +  s*  C  |-4Cos*f  j  +  | 
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where 


o  ^  A^X  ,  when  r<( 

^^io  ,  when  |fs\ 

A-  7r  f  when  ^  t  r*c  J 

Transformed  back  to  (s,  w) ,  the  solution  becomes 


where 


and 


U 


_  w*  A 

in  7T 


A  =  Von 


'  ? 


(D.7) 


zjs  3>*  (  a  -  S1 ) 

-  5  s  C  l  +  4-'SihX-^  )  +  l 


■} 


;T> .  -  ,i  >  ■  1 


v 


i  A  <-7^  when  os  w  | 

A  »  6  when  o  ^  u)  6  3.  S »  | 


A*7t  when  sJ«.coSl-|n  ^>  =  1 


Case  E: 
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To  find  u  due  to  edge  effect,  we  have  to  solve  the 
following  problem: 


By  same  transformation  formula  as  used  in  Case  B,  the 
problem  is  transformed  to 


which  has  same  solution  as  that  of 
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Placing  vortices  as  shown  in  the  figure,  we  find  that 
the  solution  will  be 


where 


—  (A,  4 

A,  ~  +3*1 


(E05) 


-»  V  43  U-  a-r  Cos  6) 
irx  -  ir  cos  ©  -1 


< a t  <  x  when  5-  <  ©  <  x  >  bs  | 


2L  <  A(  <  when  O  <  9  &  ^  5  \r  £  \ 


and 


A^-Wf — — 

(u^-o 


is  -v  nTs r  C Cos e  +  4s  ^**> e ) 

+  r  C  C05  e  -v  >T5  s 


o<  Ai<  x  when  o  <  e  6  x 

Transformed  back  to  s,w,  the  solution  of  (E„2)  will  be 


where 


U  =  - 


trot, 

nr* 


(  * 


A' 


"7T 


A,  =-  tan 


-I 


Jf  (  Cos 

i5^-2S*Cosf  ~  I 


o  <  A,  <  X  when  1  i  w  i  | 


X 


(E.6) 


and 


A. 


2£  <  a  i  <  —  7r  when  o  <  w  4  1 


J.  -l  f  4s  4-  Js  s*  (Cos,^  +>T5-siViS) 

T&h  \  - r - — ± - - 7 - -  ~~TTm 

(  -f  3 55  (  C  es  f~  4-  >TS  "sin  E*  ) 


a  <r  f\ <  x  2 


whc  rk  6  C  U)  &  “  x 


PART  III 


SHOCK  REFLECTED  FROM  SINGLE  SLOT  IN  INFINITE  SOLID  WALL 

Consider  a  shock  reflected  from  a  slot  in  an  infinite 
wall,  i.e.  from  a  strip  of  free  surface  between  solid 
walls,  which  are  infinite  in  extent,  in  both  sides0  As 
mentioned  in  the  last  part,  the  flow  characteristics 
after  reflection  at  AA'  (Fig0  3 0  1  )  will  be  three-dimen¬ 
sional  within  the  Mach  cone.  It  is  conical  first.  Then 
after  a  distance  (d  is  the  width  of  slot)  from  AA', 
the  two  waves  from  A,  A’  will  intersect  at  p  and  the  flow 
behind  this  will  start  to  be  a  superposition  of  these  two 
conical  waves.  The  flow  is  thus  no  longer  conical.  More 
farther  down-stream,  at  a  distance  md  from  AA'  the  waves 
radiating  from  A,  A’  reaches  opposite  edges  B’,  B  respect¬ 
ively.  Behind  that,  if  allowed  to  continue  unaltered, 
the  waves  would  produce  a  non-vanishing  normal  velocity 
on  solid  wall,  which  is  physically  impossible.  To  remove 
that,  one  says  two  new  waves  begin  to  radiate  from  B',  B 
with  such  a  strength  that  it  will  annul  the  non-vanishing 
normal  velocity  on  solid  wall  and  leave  the  zero  pressure 
(u  =  o)  condition  in  the  Slot  unaltered.  These  new  waves 
will  reach  c,  c*  after  another  distance  md  and  produces 
more  new  waves. 

28 
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The  flow  in  any  one  of  the  regions  in  Fig*  (3,1) 
will  be  only  influenced  by  the  waves  radiating  from  the 
points  which  are  within  the  fore  cone  of  the  particular 
region*  For  example,  the  flow  in  region  V  will  be  a 
superposition  of  the  two-dimensional  wave,  the  initial 
circular  waves  from  A,  A*  and  the  new  circular  waves  from 
B,  B',  Care  should  be  taken  in  carrying  out  the  super- 
position^  such  that  the  proper  boundary  conditions  are 
satisfied. 

Collecting  the  results  obtained  in  Part  II,  we  have: 

in  I: 

w  s  o 

V  33  o 

(3,3) 

W  a  — 

in  I  * : 


W  =  o 


60 


o  <;  a  < k  when  o  $  co  <  7c  >  x 

A  *  o  w^en  0  i  <0  *  |r  ^  mTlx  =/x 
A»-tt  when  ~  s  w  ^  a  ,  Wv'V  =x 

_  *4^  I  x-*gx  _ •  <o 


V  —  _  jrgfS  I  X->n 

I'  X  4  (W 


—  >ngx  .  to 
»mix  *  ^ 


(3.6) 


Wx  = 


TT  — 1  (Vit^K 


where 


»*-*  f 

mttx  (  \-vy\<Z>) 

^"a  ] 

( 

X  —  am  'id-'  c<^s*-i£.  3 

when 

o  £  tv  t 

m  -ax  ^  x 

when 

~  U)  —  TT 

2-  ; 

Ht  ttX  bX 

when 

O  6  u) 

»*  1 

Vntcx  -  x 

In  above 


-Or  =  J  (  ^  +  4  )*  + 

“  “  w(-*tx) 

here  o  <  U)  <- ~X 


(3.7) 


(3.8) 

(3.9) 


in  II • j 


U)  =.  <5  when  =  o  .  M  >  - 

^  2. 


co  =  x  when  =  o  ,  lJ  <  ~  ^ 


■g#  A 

rn  7f 


m-r-r)  v-  *) 


L 


5 


0 

»A 

> 

A 

* 

when 

0  <:  U)  TT, 

Olive  X 

A  =  7V 

when 

0  ^  u)  ^ 

^  / 

o>lV  at* 

A*o 

when 

I  £(j)  £7v, 

Hi®'  =x 

V  = 


4tTo( 

7T 


Co.  x 

rn  ET' 


w. 


__  4gA  X  - 

"  T  J  j«' 


Otfr 


MR/ 


<  >  ft).  _  ^  o 
•Sun  ^  7^  ^ 


(3.11) 


(3.12) 


where  J  -  W 

(  -X  -  2rt»TD"Slh  j 


when 

0  j& 

mltveX 

when 

i  & 

•34 

cO 

MRr  =*  X 

when 

O  £ 

*°*5/ 

mv“  x 

In  above 


•ok  C  4  -  4-)*  +  il 
CO  «  +ah"1  (  *4_4T") 

0<  tO-C  -7T 


(3.13) 

(3.14) 


>i  *4 

oo  *L  o  when  i  Ja^^: 

to  *7T  when  5> «  <? ,  t)  <-|- 

in  Ills 

Symbolically  the  solution  in  region  III  can  be  written 


III  «  II  +  II*  -  I 
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The  correctness  of  the  result  may  be  verified  very 
easily  by  consideration  of  the  boundary  conditions. 
For  a  point  (x,  y,  z) 


Wjl  +  M  x; 

(3.15) 

v*.  = 

+  Vx/ 

(3.16) 

W  ]jl  “ 

(3.17) 

IV*  and 

regions  further  downstream: 

To  determine  the  flow  of  the  new  wave  reflected  from 
the  opposite  edge,  a  more  complicated  problem  results, 
which  will  be  a  subject  for  further  study.  Some  discus¬ 
sions  on  this  problem  is  given  in  Part  V0 

The  pressure  coefficient  cp  along  z  axis  in  Figc 
(3.2)  is  plotted  for  regions  I  and  III  at  x  =  2  md.  See 
Fig.  5,  p.5|  and  Appendix  III. 
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PART  IV 


SHOCK  REFLECTED  FROM  MULTIPLY-SLOTTED  SOLID  WALL 

To  determine  the  flow  due  to  a  shock  reflected  from 
a  multiply-slotted  surface,  we  should  have,  in  addition 
to  what  we  did  for  free  surface  slot  in  Part  III,  a 
similar  consideration  of  superpositions  and  reflections 
of  waves  in  the  solid  wall  part.  Also  more  new  waves 
will  come  from  other  slots.  The  entire  picture  of  the 
flow  field  depends  upon  the  ratio  of  the  width  of  slot, 
d,  to  that  of  solid  surface  between  slots,  A;  and,  more-” 
over,  on  M ^  .  For  practical  interest,  that  ratio  -i 
is  relatively  small.  In  Fig.  4.1,  4.2,  ~  ,  M  &  *  JTT 

are  used. 

As  shown  in  Fig.  4.1,  4.2,  waves  radiating  from 
A,  A*  behaves  at  that  in  Part  III  first.  After  a  distance 
m/L  from  A  A’,  the  wave  from  A  (or  A')  reaches  edges 
A'  D'  (or  A  D) ,  and  if  allowed  to  continue  unaltered, 
would  produce  a  non-vanishing  u  on  the  free  surface. 

To  annul  that  u,  a  new  wave  starts  to  radiate  from  D*  (orD) 
which  will  satisfy  the  condition  w  =  o  on  solid  wall  at 
the  same  time.  A  more  complicate  flow  pattern  will 
naturally  result  from  super  positions  of  the  new  waves 
and  t  he  additions  I  waves  f  x  ojv.  neighbor!  nc  slot  f The 


situation  will  become  even  much  more  complicated  when 
we  go  further  down-stream,,  Assume  the  slots  and  the 
solid  walls  are  uniformly  spaced  with  constant  d//\.  ratio. 
Then,  as  shown  clearly  in  the  figure,  the  flow  above 
the  slot  will  be  symmetrical  about  the  vertical  perpen¬ 
dicular  plane  bisecting  the  slot  and  the  flow  above  the 
solid  wall,  symmetrical  about  the  vertical  perpendicular 
plane  bisecting  the  solid  wall  spacing.  Moreover,  we 
need  only  to  specify  the  flow  characteristic  between  the 
LL  and  RR  lines  in  Fig.  4.1.  The  flow  right  or  left  to 
that  region  is  just  a  repetition  of  this  typical  one. 

Same  as  in  Fig.  (3.1),  the  flow  in  any  region  in 
Fig.  (4.1)  will  be  a  superposition  of  waves  radiating 
from  the  points  which  are  within  the  fore  cone  of  that 
particular  region.  By  means  of  results  of  Case  A,  Part 
II,  the  flow  characteristics  of  the  regions  which  are 
not  influenced  by  waves  from  B,  B',  |S  ,  ^  or  further 
downstream  points  can  be  determined.  That  is,  we  know 
all  these  regions  which  are  above  or  forward  of  the 
extra-heavy  lines  in  Fig.  4.1 ,  4.2. 

Use  the  coordinate  system  as  indicated  in  Fig.  4.1, 
4.2.  Formulas  of  u,  v,  w  in  regions  I,  I’,  II,  II* ,  III 
are  same  as  that  in  Part  III,  i.e.  (3.3)  to  (3.17)  in¬ 
clusively,  except  that  in  (3.8),  .(3.9),  (3.13),  (3.14), 
d  should  be  replaced  by  -r. 
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In  III',  for  a  given  point  (x,  y,  z), 


^  JSi'  — 

U  -v 

4 

vxd  “°- 

Vx  + 

V2i/ 

= 

w*.  + 

In  III'^  f  III'R  and  III",  similar  formula  may  be 
written  except  that  appropriate  tsr  ,  w  should  be  used. 

The  pressure  coefficient  Cp  along  z  axis  in  Fig, 
(4,2)  has  been  plotted  for  regions  I*,  III*,  III"  at 
x  =«  ir/2  (jv  +  3d).  See  Fig.  6,  p.52.  Detail  calculation 
is  given  in  Appendix  IV.  (p.£|) 


rC'iy'i 
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PART  V 

DISCUSSION 


For  the  problem  of  a  shock  reflected  from  a  single 
slot,  there  results,  as  has  already  been  pointed  out  in 
Part  II,  pairs  of  new  waves  when  the  original  waves 
radiating  from  one  side  of  slot  reach  the  opposite  side. 
This  is  also  true  for  the  problem  of  a  flat  plate  wing 
of  very  low  aspect  ratio  (less  than  1  for  ^  ) 

at  an  angle  of  attack  ok  <>  The  pressure  distribution  of 
the  latter  problem  has  been  calculated  by  Gunn  (reference 
4)  up  to  2.5d  (d  is  the  span  of  wing)  from  the  leading 
edge,  and  by  Stewartson  (reference  5)  with  an  asymptotic 
expression.  There  arises  the  question:*  Are  these  two 
problems  similar,  and  if  so,  can  we  borrow  some  results 
of  Gunn'sor  Stewartson* s  solution  for  our  problem? 

The  differential  equation  governing  the  velocity 
potential  or  the  velocity  components  u,  v,  w  of  either 
problem  is  a  linearized  one.  Any  solution  which  satisfie 
the  differential  equation  may  be  superposed  to  form  a  ne 
solution,  provided  the  appropriate  boundary  conditions 
are  satisfied.  This  is  the  so  called  "cancellation  wim 
method.  By  this  method,  for  example,  the  lift  distribu 

*  Question  suggested  by  Prof.  W.  R.  Sears. 
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on  a  given  wing  may  be  determined  by  cancelling  excess 
lift  on  a  related  wing  with  a  known  loading;  i„e.  the 
problem  can  be  expressed  as  the  two-dimensional  wing 
problem  plus  a  cancellation  wing.  For  the  narrow  flat 
plate  problem,  the  boundary  conditions  are: 


U  *  0 

V  »  0 

W  rn  o 


(5.1) 


which  may  be  considered  as  sum  of  the  2-dimensional  flow; 


to 


JA 

m 


v»  *  o 

w,  =  -tr* 


(5.2) 


u,  =  -2sL 

*  rn 
V,  ac  O 


W|  as  -  VC<\ 


plus  the  "cancellation  wing"  problem: 
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( r.3  ) 


The  W  component  of  the  problem  is  symmetrical  about  the 
plane  of  wing  (the  x-y  plane  for  small-perturbation 
approximation ) „  Thus  we.  ha^e. 

'aw  aw  _ 

-a*  “  “  °  throughout  the  x-y  plane. 

The  problem  of  finding  W2  of  (5C3)  is  to  solve 

a.  w  _  -d  w  -i'W 

•»*'  “  °  (5.4 

with  the  following  boundary  conditions :(  Fig  5.5 ,  hex-f  page) 
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K  Wx»0  ^ 


or  just  to  solve  (5,4)  for  the  upper  half  region  of 
/c  _x  .  . 

(b„5),  due  to  °  ,  i»e„ 


Now,  for  the  slot  problem,  we  have 


which  may  be  also  considered  as  the  sum  of  2-dimensional 
flow. 


U  j  o 

v*i  *  -  (5.8) 


and  the  "cancellation  wing"  problem: 


*s  V>  “  ^ 


V*  «  o 
w>-  Vin< 


Vv  *0 

w*  -  2Vd  (5.9) 


By  irrotationality  condition,  the  problem  to  find  U2 
of  (5.9)  is  to  solve 

*  —  _  ^  2LHi 

O**  ™*  (5.10) 
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with  the  following  boundary  conditions: 


(5.11) 


compare  (5.4)  with  (5.10),  also  (5.6)  with  (5.11).  The 
differential  equations  are  identical.  The  boundary  con¬ 
ditions  would  be  also  identical  if  in  »  1  (i.e.  M  <*,  =  ) 

and  -2d  in  (5.11)  is  replaced  by  o(  .  In  other  words,  it 
seems  that  the  expression  of  W  component  of  the  narrow 
flat  plate  airfoil  may  be  used  as  that  of  u  component 
of  our  slot  problem;  except  some  changes  of  constants. 
However,  as  we  have  shown  it  clearly  in  Case  A,  Part  I, 
there  may  be  some  singularity  in  W  while  the  assumption 
of  finite  pressure  rules  out  the  possibility  of  singularity 
in  u.  The  possible  existence  of  singularity  makes  the 
solution  of  problem  not  unique.  Hence  the  u-W  corre¬ 
spondence  breaks  down  finally. 

The  above  consideration  may  be  best  illustrated  by 
the  semi-infinite  rectangular  wing  (infinite  toward  left) 
at  an  angle  of  attack  o[  with  comparison  to  Case  A  in 
Part  31.  The  W  component  of  such  a  wing  has  been  computed 


by  Gunn  (reference  4,  p.  338,  note  that  results  given  there 
are  for  wing  extended  to  infinity  toward  right)  by  means 
of  Laplace  transform.  After  changing  to  our  notation,  W  irjside. 
Mach  cone  reads  as 


+  ~  j  sih 1  ±5  +  s;h~* 


•ud  4  M  W'  3  ^ 

/  |  -  2.  c.oia-^- 

XXTot,  I—  )77^  «)_ 


-  licJSJL  c«4- 

buy  «*. 


(5,12) 


By  consideration  given  above,  the  W2  of  cancellation 
wing  is,  putting  m  =  1, 


*  ) 


(5-13> 


In  (5013),  change  <X  into  -2c(  ,  we  have  U2  of  "cancellation 
wing"  of  slot  problem  for  m=l.  Since  u_^  =  o,  the  u  com¬ 
ponent  of  slot  problem  at  M«  will  be  equal  to  U2,  i.e. 


u  _  -  ±2*  W 

Yti-K. 


(  Ccfig') 

|  \-XI7>^  Co^-^-  j 


4-jb^  coa  — 

m7V  N  ^9  (5.14) 


The  first  term  on  right  side  of  (5,14)  is  exactly 
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same  as  (A015),  p„  120  However,  in  addition,  we  have  a 
singularity  term  here 


h5 
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APPENDIX  I 


To  plot  the  constant  pressure  lines  or  iso-bars ,( Fig.  x) 
rewrite  (A ,V  t a) 


in  the  form 


to  >7  (l-hop)  Cos^r 

I-aW']  Co-S^  -!£- 


Yb^  (  l— to 9)  C.o^~- 


Give  different  value  to  cp,  i0e0  to  'fc,  we  may  plot  a 
set  of  curves  of  >*>h  vs,  Ct)„  That  results  the  Figc  2„ 
ps  48, 

To  find  the  slope  of  the  iso-bar,  introduce  } 


such  tha't 

4  Sx  =■  ^5  ) 

=  60  (  fc  ) 
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(7) 


o 

Denote  "P  =  and  by  the  formula 

3*  Co*  =  \  +  Co'S  (A)  **  |  + 

We  may  rewrite  (2)  as 

[\  U  +  -t(  »-■&  -Jv^ 71)  (8) 

or 

(j-*)U  +  *?)J  -S*('  +  r)  +  D^“S)(l4.i?)+r]  (9) 


squaring  and  collecting,  we  have 

X4  (  H-T)*  +  sM  >  +  ?)(  ^  +  iTS  “  '+ tV*T  ) 

♦  if  (  It  if )  3.  (  \*3)|  as  ©  (10) 

Hence 


^5  =  -i(  (+.rf)(f  +  3)  s^4  (  Hif)(  _ _ 


A+  C  ;  5  -  I  ,  *  •  o 

dj  o  (12) 

as  x  ( »+ 

Thus  the  slopes  of  all  iso-bars  are  zero  at  C,  except 
possibly  not  for  iso-bars  of  T  s  o  or  oo  e 
For  P  *  o,  by  (8),  wt  Ksve 

l-  HZ~TTZ  -  *  (13) 

or  3  +  3  +  <?■”*•  =  0  (1^) 

Both  (13)  and  (14)  give ~  *  o  at  C. 

For  P  =  $  ,  i0e0  cp  =  2  4/m,  we  have,  by  (8) 

WADC  IR  52-161+  c* 


which  is  parallel  to  the  boundary  surface. 
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APPENDIX  II 


To  evaluate  the  integration  of 
W  = 


to 


"Sib- 


I  —  +  2.  cot>  co  cos  -5; 


J 


c  077J*  (  I-  »9^x(  \-rn  >)*■:.) nx co) 


d  (  *>7  j 


write 


Then 


ho 


? 


J2L  =-  -5  KV)  6 


w 


4ro(  .  CO  [  Cos'fe  +  sy  6  (  a.  Cos  tp  CPS*  €  j 


7C 


^  1  "St'n’©  C  l—  slri^O  Sin  CO  ) 


-de 


-s.n  JtL 


L  (  i-siV’e  *lnxto)  (  i-Sih^e  •siVco) 


SiVo-  (  ;tCosco  Cos  -£0 


(I 


_  Cos  co  Cos  —j  T 

; - : - —  de 

sin*  to  )  _ 


2.0- oC  . 


- s.h^ 

7T  2- 


Tde  i  •»  Cos  co  j-  si'nco  P  de 

Slhx&  Sm  Ct>  J  |  +  SiVe  Siiv’’) 


l  4-  Cos  co  — 


,  CO ‘—'Si  h  <0  P  _ djS> _ 

m  <0  J  t--3rv)>’6S)ntO 


As 


t 


d  6 


S^a0- 


Co**& 


J 


do 


l  ±.s.nt«)-sih  e  JTTs^uT 


■4ah  *(4"  l  ±$l‘n  CO  4oon 
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We  have 


w  = 


^  -fah '(4 , 4 s.hU)  +an&)  -- fan (^i-si^to  +a«e) 


But 


■fan  1  t  •+  $ihw  tan&J  --fan-1  (4  i—  sinto  fan©) 


4q^' )  (-1  i+-si‘ho)  -J  '-si^Lojjane 
^  l-  i  i-  siVi'to  4av7 © 


and 


^  I  4-  Since  -  j  l  — Sli 


Sin  U) 


—  4  ( 4  >  +  ■Sin  <1)  — ,  l-S;nCO  )  n  4  *  ^  ,'r' 


60 


n  ^  \  —  Coscoj  —  2-  Sl^  2. 


Moreover 

Jr»n©  =• 


hi  *7 


>-  r>?^ 
We  have  finally 


Cot“  6 


j 


m 


w 


_  *  to 

—  -  — —  —  . 

tv  ^  K 


I-011')  _  I a  -I  \  ^  4  »>/)  (  I-  >n*?) 

n>n  7V  '  / - 

'  (  Sih  —• 
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rfl5  =  3.873 
3.85 

3.75 

3.50 

3.25 
3.00 

2.75 

2.50 

2.25 

2.00 

1.75 

4T  =  1.732 

Coefficient  is  plotted  in  Fig.  5, 


.8711 
1.622 
2.017 
2.075 
2.045 
1.975 
1.882 
1.773 
1.648 
1.508 
1.497 
P.  51. 
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APPENDIX  IV 


Refer  to  Fig.  A  „2.  Contributions  to  pressure  coef 


ficient  alone  z  axis  due  to  waves  A  and  A* 


equal 


and  contribution  due  to  waves  arid  ot  *  are  also  equal 
Thus  along  z  axis  (y  =•  o),  in  region  III' 


cp  ~  “  XT  (  ^ 


i  /  |  I  ,i  t 

XT  i  ^  ac  ^  ^  irf  ^  m  ^ 


^  f  rn 


where 

Here  A,  has  the  same  meaning  as  /\  in  (305),  except 


with  y  -  o 


w,  •  JF57  + 


CO  \  —  -la  ^  '  (  _a.  ^ 


Let  x  =  m  (  i  •— )  and  make  all  length  dimensionless 
by  dividing  with  —  as  in  Appendix  Hie  Take  =  2,50 


Then 


a,  -  +»„■'  f 

‘  l  i  -  7,** l 


Qy  i  t - ~~~~3 

t.»5  +  ** 


~  W  , 


<Jl  <  TV 


For  along  z  axis  in  region  III",  we  have„ 

F 

in  addition  to  that  given  by  (2),  the  contributions 

A 

from  waves  cK  ,  0  Denote  the  latter  by  ? 
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A  in  (3 „5j ;  except  with  y  * 

(9) 

(10) 

(11) 
(12 ) 

(13) 

Pj  _  •  pA  4  .8990 

Pp„  -  •PO  r.  3 , 1623 

f’3 ,  z*  *  46  2., 4495 

Hence  write  c  = 

P  t  m 

(0  for  z*  a  505  to  424  =?  4 

la)  for  z*  =  424  to  4To 

=  -§-A|  -4 

(  humericel  values  in  nex4  page-) 
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.#] 


J24"  =  4  0  8990 

4 

4.85 

3.1675 

4  0  75 

2.6061 

4.50 

1.9012 

4.25 

1.4927 

4.00 

1 . 2024 

3.75 

.9776 

3.50 

.7937 

JlO  =  3.1623 

.5887 

13)  for  z*  =  Jl0  to  J6~ 

3.1623 

.5887 

0 

O 

0 

co 

.5032 

.2257 

2.75 

.3845 

.3171 

2.50 

.2788 

.3523 

2.4495 

.2590 

.3555 

Coefficient  is  plotted  in  Fig.  6,  p»  52. 
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